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NONLINEAR  ANALYSIS  OF  A  RELATIVISTIC 
BEAM-PLASMA  CYCLOTRON  INSTABILITY 


I.  Introduction 

Electromagnetic  electron  cyclotron  interactions  have  been  applied  to  a 
wide  range  of  phenomena  ranging  from  radiation  and  acceleration  processes  in 
astrophysics  to  laboratory  radiation  and  acceleration  sources.  In  general  an 
electron  cyclotron  instability  results  from  a  resonance  between  the  electron's 
cyclotron  motion  with  the  Doppler-shifted  electromagnetic  wave  frequency.  One 
class  of  coherent  electron  cyclotron  instabilities  which  relies  on 
relativistic  effects  associated  with  the  electron  motion  is  called  the 
electron  cyclotron  "maser"  or  "gyrotron"  Instability.  In  this  instability 
the  electron's  cyclotron  frequency  is  a  function  of  the  electron's  energy 
(nonisochronous  rotation).  Under  certain  conditions,  i.e.,  if  the  Doppler 
shifted  radiation  frequency  is  slightly  higher  than  the  electron  cyclotron 
frequency,  the  electrons  can  experience  azimuthal  phase  bunching.  The 
resulting  electron  current  density  excites  coherent  radiation.  This  mechanism 
was  proposed  independently  in  the  late  1950's  by  a  number  of  researchers1 “**  and 
is  the  basis  of  the  electron  cyclotron  maser  radiation  source  (gyrotron)®-2® 
and  is  also  generally  believed  to  be  responsible  for  radiation  observed  from 
the  earth's  magnetic  poles, 2^"‘^1  solar  atmosphere^2-^  and  tandem  mirror 
devices. 35*36 

Under  conditions  of  resonance,  electromagnetic  cyclotron  waves  can  also 
result  in  electron  acceleration.  This  mechanism  has  been  proposed  as  a 
potential  laser  driven  acceleration  scheme®^”®^  as  well  as  a  possible 
mechanism  responsible  for  the  high  energy  electrons  observed  in  type  III  solar 
burst. 

Manuscript  appro  yd  July  11, 1985. 


In  this  paper  we  derive  a  set  of  one  dimensional,  fully  nonlinear  and 
relativistic  wave-particle  equations  which  describe  the  evolution  of  a  system 
of  particles  in  a  magnetic  field.  Both  the  beam  and  plasma  electron  dynamics 
are  treated  fully  nonlinearly  and  self-consistently  and  the  ions  are  assumed 
infinitely  massive.  This  formulation  is  used  to  study  a  beam-plasma  cyclotron 
instability  which  may  have  applications  as  a  source  of  powerful  coherent,  high 
frequency  radiation.  The  particular  cyclotron  instability  which  we  consider 
in  some  detail  is  the  high  frequency  Doppler  shifted  interaction.  This 
interaction  has  been  extensively  studied  in  the  absence  of  a  background 
plasma. 22“ 2** *3®  in  the  absence  of  a  plasma,  the  concept  which  utilizes  a  self 
resonance  condition  on  the  high  frequency  Doppler  shifted  interaction  i3 
referred  to  as  the  cyclotron-autoresonance  maser1*1  (CARM).  For  an  excellent 
review  of  this  promising  millimeter  and  submillimeter  radiation  mechanism,  see 
Ref.  19.  For  a  magnetized  relativistic  electron  beam  propagating  in  a  low 
density  plasma  the  frequency  of  this  interaction  is  u  »  Y{1q  where  Y  is  the 
relativistic  mass  factor  and  is  the  nonrelativistic  electron  cyclotron 
frequency,  and  can  substantially  exceed  the  nonrelativistic  electron  cyclotron 
frequency.  For  MeV  electron  beams,  magnetic  fields  in  the  tens  of  kG  range 
and  plasma  densities  of  -  lO1^  cm“^,  the  radiation  wavelength  is  weLl  into  the 
millimeter  regime.  For  the  Instability  to  develop  the  electron  beam  must  have 
a  nonzero  average  transverse  velocity  component  with  respect  to  the  magnetic 
field.  We  find  that  if  the  electron  beam  transverse  velocity  normalized  to 
the  speed  of  light  is  equal  to  the  inverse  of  the  relativistic  mass  factor, 


i.e.,  g  -  1/Y,  the  high  frequency  instability  can  be  optimized  in  various 
ways.  For  8  -  1/Y  an  autoresonance  exists  between  the  streaming  and  gyrating 

beam  electrons  in  a  background  plasma  and  a  propagating  electromagnetic 


wave.  The  background  plasma,  although  passive,  plays  an  important  role  in  the 
interaction  by  effectively  modifying  the  dielectric  medium  in  which  the 


electron  beam  propagates.  Also  conversion  efficiencies,  radiation  power/beam 

power,  can  exceed  20%  for  practical  device  parameters.  The  dependence  of  beam 
velocity  spread  on  the  efficiency  is  studied.  We  find,  for  example,  that  if 
f}^  -  1/Y  the  effects  of  a  finite  axial  particle  velocity  spread  has  little 
influence  on  the  efficiency,  however,  a  spread  in  the  particle's  pitch  angle 
can  significantly  reduce  the  efficiency. 

In  Section  II  the  beam-plasma  particle  orbit  equations  are  derived  for  a 
right  hand  circularly  polarized  electromagnetic  field.  The  equations 
describing  the  evolution  of  this  field  are  derived  in  Section  III.  In  Section 
IV  the  linear  regime  of  the  various  cyclotron  interactions  are  discussed  with 
emphasis  on  the  high  frequency  mode.  Sections  7  and  VI  discuss  the  effects  of 
a  beam  velocity  spread  on  the  instability  and  in  Section  VII  an  analytical 
expression  for  efficiency  is  obtained.  The  results  of  numerically  solving  the 
system  particle-wave  equations  are  presented  in  Section  VIII  for  laboratory 
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II.  Particle  Dynamics 

In  this  section  the  relativistic,  nonlinear  particle  orbit  equations  are 
obtained  in  a  convenient  representation.  These  particle  orbits  are  used  in 
Section  III  to  derive  the  appropriate  driving  currents  for  the  electromagnetic 
field.  The  current  densities  are,  in  turn,  used  to  derive  a  system  of  self- 
consistent,  fully  relativistic  equations  describing  the  temporal  evolution  of 
the  wave. 

The  transverse  vector  potential  associated  with  the  right  hand  circularly 
polarized  electromagnetic  field  is  represented  by 


I 


P 


A(z, t)  -  A(t)  (sin  <fr(z,t)e  ♦  cos  $(z,t)e  ),  (1) 

**  ^  y 

where  A(t)  is  the  temporally  varying  amplitude  of  the  circularly  polarized 
wave,  <p(z,t)  -  kzz  -  u(t')  dt*  is  the  phase  of  A(z,t),  kz  is  the  wave 

number,  and  u(t)  is  the  temporally  varying  frequency. 

The  relativistic  orbit  equations,  for  electrons  in  the  fields  given  by 

a 

Eq.  (1)  together  with  the  axial  magnetic  field  B  -  Be,  are 

**“  o  z 

B  ■  (Q0  x  "  Ie|  ( "C  1  3A/3t  ♦  g  x  (VxA)/YmQ)  (2) 


where  p  «  Km  v  is  the  particle  momentum  vector,  Q.  -  a_e_,  0  *  iei  B  /m  c  is 

O  “*  y  v  4  Q  1  1  O  O 

2  1  /2 

the  nonrelativistic  electron  cyclotron  frequency,  Y  -  (1  *  (g/mQc)  ) 
is  the  Lorentz  factor,  mQ  is  the  electron's  rest  mass,  e  is  the  electronic 
charge  and  the  dot  denotes  a  total  time  derivative.  Substituting  (1)  into  (2), 
the  three  components  of  the  momentum  equation  become 


Px  ■  *  Py  ■  C(<u(t)  -  v„k„)  cos  $(z,t)  -  f(t)  sin  *(z,t)], 


z  z 


(3a) 


(3b) 


py  -  “  Px  ♦  —  C(«(t)  -  vz kz)  sin  $(z,t)  +  r(t)  cos  <p(z,t)], 


Pz  -  -  | e |  kz  [px  cos  *(z,t)  -  p  sin  <p(z,t)], 


where  r(t)  «  A  (t)  3A(t)/3t  is  the  nonlinear  growth  or  decay  rate  of  the 
electromagnetic  field  and  vz  -  pz/YmQ.  A  convenient  representation  for  the 
transverse  particle  momenta  is  p  -  p  cos  0,  and  p  -  p  sin  8,  where 

x  1  y  i 

p^  -  p^(zQ,  9q»  E0»  t)  is  the  magnitude  and  9  -  8(zq,  ®0  E0»  is  the  Phase 

of  the  transverse  particle  momentum  at  time  t  having  the  initial  values  zQ,  9Q 
and  qq.  Using  the  transverse  momentum  representation,  the  orbit  equations  in 
(3)  become 


P  -  -  l-l---  C(ui  -  p  k  /Ym  )  cos  ip  -  r  sin  ip], 
JL  ®  z  z  o 


ip  -  ~(to  -  pk/Ym  -  Q  n) 
z  z  o  o 


+  ~  C(«  “  p_k  /Ym  )  sin  ♦  ♦  r  cos  ip], 

c  p  z  z  o 

1 


lel  Akz 

ITT-  p,  003 

O  J- 


where  ip  -  e(t)  +  <p(z,t)  i3  the  relative  phase  between  the  particle  transverse 
momentum  and  the  electromagnetic  wave.  The  rate  of  change  of  Y  in  this 
representation  is 

|  ®  |  Ap 

Y  - — — (u  cos  ip  -  f  sin  ip).  (5) 


Ill .  Self-Consistent  Evolution  of  Fields 


The  one  dimensional  wave  equations  for  the  vector  potential  A(z,t)  is 


-  1.11.)  A 

c2  3t2  ' 


(6) 


where  J(z.t)  is  the  response  particle  current  density  which  will  be  determined 
by  taking  appropriate  averages  over  the  particle  orbits.  Substituting  (1)  into 
(6),  the  wave  equation  can  be  cast  into  the  form 


a(t)  cos  <fr(z,t)  ♦  B(t)  sin  $(z,t)  -  -4*3  Jx(z,t),  (7a) 

-a(t)  sin  $(z,t)  ♦  6(t)  cos  <fr(z,t)  -  -4irc  Jy(z,t),  (7b) 

where  a(t)  -  (2ru  +  3w/3t)A  and  B(t)  *  («2  -  c2k2  -  r2  -  3r/3t)A. 

z 

To  obtain  a  and  8  explicitly  in  terms  of  the  response  current  density  we 
multiply  (7a)  and  (7b)  by  cos  <p  and  sin  $.  Upon  integrating  the  resulting 
equations  over  z  from  0  to  2*/k_  we  obtain 

Zi 

r2ir/kz 

o(t)  »  "2kz°  J0  dz(Jx(z,t)  cos  $(z  ,t)  -  Jy(z,t)  sin  <f(z,t)),  (8a) 
r2ir/kz 

8(t)  ■  “2kzc  J0  dz(Jx(z,t)  sin  ^(z,t)  +  Jy (z ,t)  cos  <j»(z,t)),  (8b) 

In  Eqs.  (8),  the  spatial  integral  is  taken  from  0  to  2ir/kz  rather  than  from  0 

to  <».  The  reason  for  this  is  that  we  are  considering  solutions  for  the  fields 
that  are  spatial  periodic  with  a  period  equal  to  a  wavelength  2ir/kz  and  with 
temporally  varying  amplitudes  and  frequencies.  The  temporal  evolution  of  a  and 
8  completely  characterizes  the  electromagnetic  wave.  To  obtain  the  temporal 


dependence  of  a,  and  B,  and  hence  the  wave  amplitude  and  frequency,  the  self- 
consistent  current  densities  are  needed. 

The  particle  current  density  can  be  expressed  in  terms  of  the  particle 
velocities  averaged  over  the  initial  particle  distribution  function.  The 
average  current  density  can  be  shown  to  be  given  by 

J(z,t)  -  -|e,nb  I  d3po  /  dzQ  f0(z0,E0.  tm°)  Y<V  V 


5(z-z(zo,  t)), 


(9) 


where  nK  is  the  average  ambient  particle  density,  z(z  ,  p  ,  t)  and 

u  O  *"0 

v(zo,  bq,  t)  -  denotes  the  axial  position  and  velocity  of  the  particle 

with  initial  position  zQ  and  initial  momentum  vector  pQ,  fQ(z0,  p0,  t-0)  is  the 

initial  Vlasov  distribution  function  and  the  integrals  are  performed  over  all 

initial  momenta  and  positions.  The  initial  distribution  function  fQ  is  assumed 

to  be  spatially  symmetric  about  the  z  axis  and  normalized  such  that 

|  fQ(zo,  bo>  t-0)  d3pQ  -  nb(zQ,  t-0),  where  nb(zQ,  t-0)  is  the  initial  beam 

particle  density  and  dp  -  p  dp  dp  d0  .  In  obtaining  the  current 

o  oj_  oj_  oz  o 

density  representation  in  (9)  the  actual  particle  distribution  function  was 
written  as  a  sum  of  the  individual  particle  distribution  function.  For  a  large 

number  of  particles  the  sum  can  be  replaced  by  integrals  over  the  initial 

momentum  and  position  of  the  particles.  Generalizing  the  response  current  in 
(9)  to  more  than  one  particle  species  evolves  summing  (9)  over  the  various 
particle  distributions.  Substituting  (9)  into  Eqs.  (8)  together  with  the 
representation  for  the  transverse  momentum,  and  performing  the  specified 
integration  over  z  yields 


7 


/a(t) \ 
6(t) 


4it  |  e  |  n^c 
m 


[2it  //  p  dp  dp  g  (p  ,p  ) 
Oj_  °j_  02  o  Oj.  02 


(cos  r|» 
sin  ip 


],  (10) 


where  ip  -  ip(z  ,  p  ,  t),  the  average  is  defined  as  <•••>.  -  4—  f2ir  dip  , 

o  co  ^  2l  '0  o 

ip„  -  <p„  +  9„  and  g  (p  ,  p  )  -  f  (z  ,  p  ,  t-0)  is  the  initial  distribution 
o  to  o  0  0j_  02  o  o  fco 

function  which  is  assumed  to  be  axially  symmetric  and  spatially  homogeneous. 
The  initial  Vlasov  distribution  is  normalized  such  that 

2n  fj  p0^dp0^dpozg0(p0^*Poz)  “  1  •  In  obtaining  (10)  we  used  that  relationship 

2ir/k_  J  f2ir 

kz  ^  dzQ  -  JQ  dipQ.  In  general  the  integration  over  zQ  in  Eqs.  (10) 

should  include  all  initial  particle  positions.  However,  in  Eqs.  (8)  the 
integration  over  z  extends  from  0  to  2ir/kz  and  since  the  system  of  particles 
and  fields  are  periodic  in  z,  we  may  take  the  zQ  integration  range  in  (10)  to 
be  from  0  to  2ir/kz. 


Normalized  Particle-Wave  Equations 

The  expressions  in  (10)  together  with  the  orbit  equations  in  (^ )  describe 
the  fully  nonlinear  evolution  of  the  fields  represented  by  (1).  Introducing 
normalized  variables,  the  orbit  equations  in  (M)  become 


3u 

—1 

3t 


3ip 
3t  " 


3u 

3t 


-  u  k 

-  -a[(io - - — )  cos  ip  -  r  sin  ip], 

_  _  u  k  _ 

-Afl  +  —  [fu - - — )  sin  -p  +  r  cos  ip], 

1 

u 

-  -a  k  cos  ip, 

z  Y 


(11a) 


(lib) 


(11c) 


and  the  field  equations  in  (10)  become 
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u 

p  dp  dp  g  <  cos  ip  >  1,  (12a) 

Oi  Ol  02  O  Y  \|» 

o 

_2  u 
a  —  a>  [  2tt  /  f  p  dp  dp  g  <  -^  sin  \l>  >  ],  (12b) 

b  oj_  o_l  oz  o  Y 

where  u  -  p  /me,  u_  -  p_/m.c,  t  -  £2  t/<Y  >,  <Y_>  is  the  average  initial 

L  L  O  £»  £»  \j  oo  y 

value  of  Y,  a  -  | e |  A/m0c2,  u  =■  ui/(no/<YQ>) ,  kz  ■  ck2/(QQ/<Yo>) , 

k  -  ck  /(Q  /<Y  >),  F  -  r/(Q  /<Y  >),  Afi  -  w  -  u  k  /Y  -  <Y  >/Y, 
J_J_00  00  zz  o 

—  2 

wb  *  (^/(q  /<Y  >)  and  “u  *  ( | e i  nb/mQ)  is  the  plasma  frequency.  For 
completeness  we  include  the  normalized  form  of  (7) 

u 

dY  |  -  - 

—  *  a(u  cos  i|>  -  r  sin  i|>).  (13) 

dt  Y 

The  initial  conditions  for  Eqs.  (11),  (12)  and  (13)  are  u  (t-0)  ■  u  , 

■  1  °1 

uz(t-0)  ■  uqz,  Y(t-O)  -  Yq,  ip ( t  —0 )  -  i|>o,  w(t»0)  -  ti>o  and 
F(t-O)  -  0. 

Up  to  this  point  we  have  considered  a  general  distribution  of  particles 

interacting  nonlinearly  with  electromagnetic  and  electrostatic  waves.  The 

present  formulation  can  treat  several  particle  species  such  as  a  beam  and 

background  plasma  by  simply  taking  the  initial  Vlasov  distribution  function 

g  (p  ,P  _)  to  be  a  sum  over  all  the  particle  species.  In  this  approach  all 
O  0 1  oz 

the  particle  species  are  treated  fully  nonlinearly  and  self -consistently. 
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Linear  Growth  Rates  for  Beam-Plasma  Cyclotron  Interactions 


The  linear  electromagnetic  dispersion  relation,  in  unnormalized  units, 
for  a  cold  magnetized  electron  beam  propagating  and  gyrating  within  a  cold 
stationary  plasma  can  be  derived  from  the  fully  nonlinear  equations  in  (11) 
and  (12).  The  details  of  the  derivation  will  not  be  given  here,  but  can  be 
found  in  Ref.  (3*0.  The  electromagnetic  dispersion  relation  is 
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(14) 


(to  -  V  k  -  Q  /Y  ) 

oz  z  o  o 


The  left-hand  side  of  (14)  represents  the  electromagnetic  dispersion  relation 
for  right-hand  circularly  polarized  waves  in  a  cold,  stationary  plasma  with 


plasma  frequency  cop.  The  right-hand  side  of  (14)  represents  the  coupling  term 


due  to  the  propagating  and  gyrating  beam.  The  beam  is  taken  to  have  a  plasma 


frequency  tow,  axial  velocity  v„  and  transverse  velocity  v  -  cB  . 

u  oj_  oj_ 

In  general  there  are  three  unstable  interaction  frequencies  associated 


with  the  dispersion  relation.  Figure  1  shows  the  general  features  of  the 
dispersion  relation  and  the  three  possible  intersection  frequencies  for  waves 
propagating  along  the  magnetic  field.  In  this  section  we  will  be  primarily 
concerned  with  finding  the  linear  growth  rates  for  the  medium  and  high 
frequency  modes,  labeled  and  u2,  respectively.  To  solve  for  the  growth 


rates  we  set  w  -  u  ♦  6co,  where  w  -  v  k  +  n  /Y  and  is  the  Doppler 

o  oozzoo 


shifted  cyclotron  frequency  and  |wo|  >>  |6u|.  The  dispersion  relation  reduces 


(15) 
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where  we  have  set  o>  -  vz kz  equal  to  Q0/Yq  on  the  right-hand  side  of  (14)  and 
have  taken  the  Doppler  shifted  cyclotron  frequency  w0  to  lie  on  the 
electromagnetic  dispersion  curve,  hence  wQ  satisfies  to2  -  c2k|  -  (^^/(wq  - 

V  -  °* 


Case  1  (medium  frequency  interaction,  w  -  u1  >  Qq) 

In  the  absence  of  the  beam,  -  0,  the  frequency  of  the  electromagnetic 
wave  in  the  background  magnetized  plasma,  for  waves  slightly  above  the  cut-off 
frequency,  is  given  by 


w  -  w 
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(16) 
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where  w _ »  Q_  (1  ♦  <o/Q  )  is  the  cut-off  frequency,  ick  |  <  u  and  we  have 

co  o  p  o  1  z 1  co 

assumed  u  <<  Q  . 

P  o 

In  this  limit  the  dispersion  relation  in  (15)  reduces  to 


a3  u2  a  n  -  a2  a  n2  a2 

( - 2 — _)  5u  .  _S  ;  °  ^  °1  °  -  -|i  JL], 

ft  6ai  +  id  O  6u> 
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where  -  <u-u>0  and  ojq  is  given  by  (16).  Assuming  that  |6w|  <<  Yo8^ftQ/2  we 


find  that 
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(18) 


(m  +  fl  5m) . 
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From  (18)  we  find  that  the  growth  rate  for  a  wave  propagating  along  the 


magnetic  field  with  a  frequency  near  the  cut-off  frequency, 


m  -  Q  (1  ♦  m2/fl2)  *  Q  ,  is 
co  o  p  o  o 
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for  r  <<  m  /G  or 
P  o 


for  r  »  m2/Q  . 
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This  result  should  be  compared  to  the  usual  growth  rate  associated  with  the 
cyclotron  maser  instability.  That  is,  taking  uip  -  0,  kz  -  0  and  introducing  a 
transverse  wave  number  k^,  it  can  be  shown  that  the  usual  cyclotron  maser 
growth  rate  is 


2 


Comparing  (19a)  with  (20)  we  note  that  the  effect  of  the  background  plasma  is 

equivalent  to  an  effective  transverse  wave  number.  The  effective  transverse 

wave  number  introduced  by  the  background  plasma  is  not,  however,  i»>co/c,  as 

would  be  expected,  but  is  substantially  smaller,  due  to  the  presence  of  the 

2  2  2 

nearby  cyclotron  mode,  and  equal  to  (w  /u>  )ar  /c  -  m  /(fi^c). 
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Case  2  (high  frequency  interaction,  u  -  w2  »  ao^ 

22  21/2 

For  the  case  where  <o9  >>  we  find  that  u  -  (c  k  ♦  «  ) 

^  o  o  z  p 

2  2  2 

*  ck  (1  +  «  /(2c  k  )),  where  u.  »  >>  Hence  the  dispersion  relation  in 

2  p  Z  o  w  p 


(15)  reduces  to 
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Since  a)  -  v  k  +  Q  /Y  -  (ok  ♦  w*) 
o  z  z  00  z  p 

interaction  is  given  by 


the  wave  number  for  this 


ck  -  y2  -d  +  H  "  “nYo/CQoYL))1/2/8oz) 
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where  Y  -  (1  -  82  ) 1/2  and  8  -  v  /c.  Using  (22)  we  find  that  the 

oz  oz  oz  oz 

interaction  frequency  o>0  is 


o  r.  ^  „2  .2  /fl  ^  ,,  2v,2/._2„2  .1/2., 

oi  -  —  C1  +  Y_8^(B„  *  0  “  «_Y  /(Q  Y  )  )]. 
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For  a  highly  relativistic  electron  beam  in  a  tenuous  plasma,  i.e.,  YQZ  »  1 

and  u  Y  /(Q  Y  )  «  1,  the  interaction  frequency  is  2Y  2  times  higher  than  the 
p  o  o  oz 

2  2  2 

relativistic  cyclotron  frequency,  that  is,  u>  -  2Y  0  /Y  -  2Y  G/O+Y  B_  ). 

J  ’0  ozoo  00  o  oj_ 

Two  particularly  simple  limiting  cases  for  the  high  frequency  interaction  can 
be  considered. 


Case  2a  (i5on  >>  ajfe2  /2|Q_/Y_  -  8  2  <u_i) 

—  .  1  '  p  0_L  1  o  O  Oj  o' 

In  this  limit  the  dispersion  relation,  for  a  highly  relativistic  beam  in 

2  2  2 

a  tenuous  plasma,  reduces  to  6u  -  (1  -  B  -  6„  )/2Y^  where 

r  b  oz  oj^  o 
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w0  -  2YQ2  Q0/Yq.  This  high  frequency  interaction  is  unstable 
2 

if  8Q^  +  8qz  >  1  and  the  growth  rate  is  given  by 


/W  °i  02 


Note  that  (24)  is  valid  for  wave  frequencies  satisfying  (23).  no  frequency 
mismatch.  For  a  more  complete  description  of  the  linear  theory  in  thi3  limit 
see  Ref.  (15). 


Case  2b  (i5wi  <<  w?82  /2|Q  /Y_  -  82  u. i ) 
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In  this  limit  the  dispersion  relation  in  (21)  reduces  to 
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The  linear  growth  rate  associated  with  this  high  frequency  mode  is 


2  2 
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r  2  l4u  y  °oi 


In  this  section  the  linear  regime  associated  with  the  beam-plasma  cyclotron 
interaction  has  been  briefly  discussed  for  the  medium  and  high  frequency 
interaction.  The  remainder  of  this  article  will  deal  primarily  with  the 
effects  of  beam  velocity  spreads  and  saturation  efficiency  of  the  high 
frequency  interaction. 


V.  Effects  of  Velocity  Spreads 


A  somewhat  simplified  but  useful  way  of  considering  the  role  played  by  a 
beam  energy  spread  is  to  examine  the  cyclotron  resonance  denominator  term, 

Q  -  a)  -  vQZkz  -  8Q/Y0  in  the  linear  dispersion  relation.  For  a  beam  with  zero 
energy  spread  in  both  the  longitudinal  and  transverse  directions  the  resonant 
denominator  in  the  linear  interaction  regime  is  nonzero  and  equal  to  the 
linear  frequency  shift  Aw.  ,  that  is  8(v  ,v  )  -  «  -  v  k  -  8  / Y  «  Aw.  , 

L  OZ  Oj_  OZ  Z  O  O  L 

where  [Acu^i  <<  u.  If  the  beam  electrons  have  a  spread  in  longitudinal  and 

transverse  velocities  denoted  by  Av_  and  Av  respectively,  the  cyclotron 

2  1 

resonance  term  will  in  turn  have  a  spread,  A8.  If  |A8|  is  small  compared  to 
the  magnitude  of  the  linear  frequency  shift,  |Aa^|,  the  interaction  will 
evolve  as  if  Avz  and  Av^  were  zero.  If  on  the  other  hand  |AQ|  is  large 
compared  to  iAii^i,  the  effect  on  the  instability  of  the  finite  Avz  and  Av^ 
will  be  significant.  These  simple  physical  arguments  can  be  quantified  by 
expanding  8(v  .v  )  about  v  ,  v  ,  where  v_  -  v„  ♦  Av_  and 

Z  ^  OZ  Q  |  At  s it 

v  -  v  ♦  Av  .  Expanding  8  gives  8(v  ,v  )  -  8(v  ,v  )  +  A8, 

X  X  z  oz  oj_ 

where  A8  -  (38/3v  )Av  +  (3a/3v  )  Av  ,  38/3v  -  -k  ♦  (8  /Y2)3y/3v  - 

zz  zzo  z 

-k  +  Y  8  v  /c2  and  38/ 3v  -  (8  /Y2)  3Y/3v  -  Y8  v  /c2.  The  spread  in  the 

z  oz  X®  X®X 

cyclotron  resonant  term  due  to  velocity  spread  is  therefore 


A8  -  (-ck  M  a  J  ♦  Y  8  8  a)  Av  /c,  where  a  -  Av  /Av  .  The  effects  of 
z  o  o  oz  o  o  oj_  z  '  X  z 

velocity  spread  on  the  evolution  of  the  instability  can  be  neglected  if 


lckz  "  W8oz  +  *o±  a)l*Vc  <<  I  AuL I  * 


The  effects  of  a  longitudinal  velocity  spread  can  be  removed  to  lowest  order 

if  Y  8  B  -  ck  .  For  Y  8  S  -  ck  condition  (27)  becomes 
o  o  oz  z  o  o  oz  z 


The  spread  In  the  relativistic  mass  factor,  Y,  associated  with  the  velocity 

spread  is  AY  -  Y-Y  -  Y3(8  Av  +3  Av  )/c. 

o  o  oz  z  o_l  J. 

For  a  particular  choice  of  parameters  an  interesting  possibility  arises 

for  significantly  relaxing  the  effects  of  a  longitudinal  beam  velocity  spread 

on  the  instability.  If  we  consider  the  high  frequency  interaction  regime 

(Case  2)  for  a  highly  relativistic  electron  beam  in  a  tenuous  plasma  the 

2 

interaction  frequency  is  given  by  u  •  ck  »  2Y  Q  /Y  .  Using  (23)  and  the 

O  2  OZOO 

condition  Y  Q  3  -  ck  *  u  we  find  the  necessary  value  of  3  required  to 

o  o  02  2  o 

remove  the  effects  of  a  longitudinal  velocity  spread  on  the  high  frequency 


mode  is  8  -  1/Y  .  Taking  8 

oj_  o  oj_ 

condition  in  (28)  becomes 


1/Y  ,  the  transverse  velocity  spread 


(29) 


Note  that  for  3 

°1 


1/Y  the  interaction  frequency  is  uQ 


YqQ0  and  the 


condition  for  Case  2b  (see  Eq.  (25))  is  automatically  satisfied  since 


2 

Q  /Y  -  8  m  vanishes, 
o  o  ox  o 


VI.  Beam  Quality  and  Velocity  Spread 


The  measure  of  quality  of  a  spatially  uniform  electron  beam  propagating 


along  and  gyrating  about  a  uniform  magnetic  field  is  given  by  two  independent 


quantities.  In  this  paper  we  will  choose  these  two  independent  quantities  to 


be  the  fractional  half-width  spread  in  the  relativistic  mass  factor  V  and  the 


velocity  pitch  angle  6.  where  8-8/8  .  That  is,  the  measure  of  beam 

°  O  Oj_  oz 


quality  is  determined  by  A Y/yq  and  60/0q  where  AY/Yq  and  A6/8q  are  positive  and 


assumed  small  compared  to  unity.  Figure  2  shows  the  velocity  distribution  for  a 


beam  having  a  50  and  AY  spread.  In  Fig.  2  the  beam  electrons  are  assumed  to 


uniformly  fill  the  crossed  region.  The  full  spread  in  8Z  and  8^,  about  the 


average  values  8rt_  and  8  ,  is  respectively  2A8_  and  2A6  .  Using  the  symbol 


oz  -  -0j/  r -  ' - z - j/  ° - ' - 

6  tc  denote  the  variation  in  a  quantity  from  its  average  value,  we  find  for 


individual  electrons  that  58  -  9  58  ♦  58  -  8  (5Y/Y  +  Y250/8  )/Y2  and 

j_  o  o  o  o  o  o 

58,  -  58  -  9  50  -  (6Y/Y  -  Y20260/0  )/Y2  where  58  -  5Y/Y3  is  the  variation  in 


o  o  o  o 


the  total  electron's  velocity  and  Yq  >>  1.  Mote  that  the  variations  58^, 


68  ,  68,  50  and  6Y  can  be  of  either  sign.  For  0  -  1/Y  and  Yrt  »  1  the 

z  oj^  o  o 


average  pitch  angle  is  0  -  1/Y  and  the  beam's  half-width  spread  in  8  and  8 

o  o  ^  z 


is  given  by 
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The  expressions  in  (30)  give  the  half-width  spread  in  the  transverse  and 


longitudinal  velocities  in  terms  of  AY/Y  <<  1  and  A0/0  <<  1  for  Y.  >>  1. 

o  o  o 


Combining  (30)  and  (29)  we  find  that  the  fractional  half-width  spreads  AY/Y 


and  A0/0O  can  be  neglected  in  the  high  frequency  radiation  excitation  process 
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Substituting  the  linear  frequency  shift  Acu^  -  (Wpi^S^/C^Y^ )  i.e., 

the  real  part  of  5u>  in  the  solution  of  (25),  into  (31)  gives  for  the  quality 
requirements  on  the  electron  beam 
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(32) 


where  we  have  set  0  -  1/Y  and  u  -  Y  Q  .  The  inequality  in  (31)  shows  that 

oj.  o  o  o  o 

a  beam  pitch  angle  spread  is  more  detrimental  to  the  wave-particle  interaction 
than  a  beam  energy  spread  for  relativistic  electron  beams.  This  point  will  be 
demonstrated  in  Section  VIII. 
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VII.  Efficiency  and  Saturation  of  the  High  Frequency  Interaction 

A  similar  line  of  reasoning  as  used  to  determine  estimates  of  velocity 
spread  effects  can  be  used  to  estimate  saturation  efficiencies.  If  electron 
trapping  is  responsible  for  saturating  the  instability  one  expects,  for  a  beam 
without  an  initial  velocity  spread,  that  at  saturation  the  resonance  term 


becomes  Q(v  ,v  )  •  w  -  v  k  -  Q  /Y  ■  -f  iw.  ,  where  v0_  and  Y_  is  the  axial 
SZ  3j_  SZ  Z  0  S  L.  sz  s 

beam  velocity  and  relativistic  mass  factor  at  saturation  and  0  <  f  <  1  is  a 
phenomenological  factor  related  to  the  fraction  of  trapped  beam  electrons  in 


the  cyclotron  wave. 

We  now  expand  Q(v  ,v  )  about  its  initial  value  giving 

SZ 


where 
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and  Q(v  ,v  )  -  u  -  v  k  -  Q  /Y  -  Aw.  .  In  obtaining  (34),  the  relation 

OZ  Oj_  OZ  Z  O  O  b 

dY/dpz  -  (a/k2)m0c2  was  used,  this  relationship  can  be  obtained  from  the  orbit 
equations  In  (4)  and  (5).  From  (33),  together  with  J2(v  ,v  )  -  Aw,  ,  we  find 

OZ  OJ^  L 

2  2 

that  Yq  -  Yg  -  Yq(1  ♦  f)  AwLnQ/Wp  where  Aw^  is  the  linear  frequency  shift 

given  by  the  Re(6w)  in  (25).  Defining  the  saturation  efficiency  as 

n  -  (Y  -Y  )/ ( Y  -1),  we  find  for  Y0  »  1  that 
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The  approximate  efficiency  expression  in  (35)  is  appropriate  for  the  case 


where  6Q^  *  1/YQ,  i.e.,  and  therefore  is  invalid  for  very  large  and  very  small 

values  of  up  as  stated  in  Section  IV  under  Case  2.  The  efficiency  expression 

1  /3 

in  (35)  scales  as  1/Y  or  as  the  cube  root  of  the  radiation  wavelength.  A 
comparison  of  (35)  with  the  numerically  computed  efficiencies  is  given  in 
Section  VIII. 
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VIII.  Numerical  Illustrations 

In  this  article  we  are  interested  In  the  possibility  of  applying  the  high 
frequency  interaction  process  to  the  generation  of  high  power  millimeter  and 
submillimeter  radiation.  To  this  end  the  coupled  wave  particle  equations  in 
(11)  and  (12)  are  integrated  numerically  for  the  high  frequency  interaction 
which  occurs  at  u  *  *  TyPically  a  few  hundred  particles  per 

wavelength  are  integrated  for  both  the  beam  and  plasma  particles.  Note  that 
in  our  formulation  of  the  problem  only  the  dynamics  of  those  particles  in  a 
single  wavelength  are  needed. 

Figure  3  shows  the  efficiency  and  growth  rate  as  a  function  of  normalized 
background  plasma  frequency  for  a  1 .5  MeV  (y  =  4)  and  3.5  MeV  (Y  -  3) 
electron  beam.  The  transverse  beam  electron  velocity  is  BQ^  -  1/YQ  -  0.25, 
0.125  and  the  beam  plasma  frequency  is  uip  -  0.025  Qq.  For  an  external 
magnetic  field  of  BQ  -  40  kG,  the  beam  current  flux  i3  approximately  500 
A/cm  .  The  interaction  frequency  Is  at  -  YQao,  corresponding  to  a  wavelength 
of  X  -  0.67  mm  for  Yq  -  4  and  A  -  0.3^  mm  for  Yq  -  8.  The  efficiency  in  Fig. 

3  reaches  a  maximum  value  of  -  1 8%  at  a  plasma  frequency  of  oip  -  0.2  Qq 


(nr 


p  -  2  x  10^  cm  3).  The  calculated  saturation  efficiency  given  in  Eq.  (35) 

for  these  parameters  is  found  to  be  n  *  1 4  ( 1  f )  %  where  f  is  between  0  and  1,  and 

is  in  reasonable  agreement  with  the  numerical  efficiency  of  18%.  The  growth 

-3  -2  -1 

rate,  at  maximum  efficiency  for  •  4  is  r  •  2  x  10  SWc  «  4.7  x  10  cm 

which  corresponds  to  an  e-folding  length  of  Lg  -  1/r  *  22  cm. 

Figure  4  shows  the  effect  on  efficiency  for  a  beam  with  energy  spread, 

AY/Y  and  pitch  angle  spread  A0/0Q.  The  other  parameters  for  this  figure  are 

8  ■  1/Y  ,  ■  0.2  Q  and  m.  -  0.025  0  . 

oj_  o  p  ob  o 


v-  v  •  ■  f.  .\  •  V  -  *’•  •*.<■* 
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In  Fig.  5  the  saturation  efficiency  iT  shown  as  a  function  of  beam 

energy.  In  this  figure  the  transverse  beam  electron  velocity  is  -  1/Yq 

while  the  beam  density  (current  flux)  and  background  plasma  density  are  the 

same  as  for  Fig.  4.  A  noteworthy  feature  in  the  figure  is  the  relatively  slow 

falloff  in  efficiency  with  increasing  beam  energy.  The  wavelength  of  the 

excited  wave  is  0.67  mm  for  Y  -  4  and  84  pm  for  Y  -  32. 

o  o 

IX.  Discussion 

We  have  analyzed  a  beam  plasma  cyclotron  wave  instability  which  may  have 
application  as  a  high  efficiency,  high  frequency  radiation  source.  Although 
our  analysis  is  somewhat  general  we  have  concentrated  on  the  high  frequency 
Doppler-shifted  interaction  since  efficient,  high  frequency  sources  are  more 
challenging  to  develop.  The  nonlinear  coupled  wave  particle  equations  were 
numerical  solved  and  various  plots  of  the  radiation  efficiency  as  a  function 
of  background  plasma  frequency,  beam  energy  spread  and  beam  particle  pitch 
angle  spread  are  presented.  The  conversion  efficiency  in  the  high  frequency 
plasma  maser  source  could  be  enhanced  by  appropriately  contouring  the  external 

1  Q 

magnetic  field  as  a  function  of  axial  position.  ’  Spatially  contouring  the 
external  magnetic  field  in  order  to  enhance  efficiency  has  been  successfully 
employed  in  conventional  maser  (gyrotron)  devices.1^ 

By  choosing  the  average  transverse  beam  velocity,  vq^,  equal  to  c/YQ  the 
deterimer.tal  effects  of  a  spread  in  the  axial  beam  velocity  Av_  can,  to  lowest 
order,  be  removed.  We  have  shown  that  high  efficiencies,  tens  of  percent,  can 
be  achieved  even  for  electron  beams  with  large  AY/YQ  spreads  and  modest  A6/eo 
spreads,  see  Fig.  4.  Both  analytical  and  numerical  results  show  that  the 
saturation  efficiency  falls  off  gradually  as  the  beam  energy  is  increased,  see 
Fig.  5. 
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Although  not  contained  in  the  present  paper  our  analysis  shows  that  the 


present  radiation  mechanism  may  be  more  applicable  to  an  oscillator 
configuration.  The  high  frequency  Doppler-shifted  mode  the  background  plasma 
merely  serves  to  modify  the  index  of  refraction  of  the  electromagnetic  wave, 
making  it  less  than  unity.  In  principle,  the  background  plasma  can  be  removed 
and  replaced  by  a  propagating  wave  having  a  transverse  spatial  gradient  such 
as  waves  in  an  open  resonator  configuration.  An  open  resonator  configuration 
has  a  number  of  important  characteristics  and  advantages  over  a  closed 
resonator  or  waveguide  configuration. 

In  an  open  resonator  transverse  mode  selection  is  a  natural  consequence 
of  the  larger  diffraction  angles  associated  with  higher  order  transverse 
modes.  Higher  order  transverse  modes  tend  to  propagate  out  of  the  resonator 
and  thus  have  lower  Qs.  Longitudinal  mode  selection  can  be  achieved  by  using 
a  Fabry-Perot  interferometer.  Furthermore,  since  the  surface  area  of  the 
mirrors  can  be  made  large,  open  resonators  can  operate  at  extremely  high  power 
levels. 

A  possible  configuration  for  a  plasma  cyclotron  maser  oscillator 
utilizing  an  open  resonator  cavity  is  shown  in  Fig.  6.  Radiation  which  builds 
up  within  the  resonator  is  extracted  through  a  small  opening  in  the  right 
mirror. 
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Figure  6  Shows  a  possible  configuration  for  a  plasma  cyclotron  maser 
oscillator  source. 
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